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I. INTRODUCTION 



Cosmological perturbation theory 0, 0] has become the standard tool to study inflation and its obseryational 
consequences for the Cosmic Microwave Background (CMB) and the formation of large scale structure 3J. Lin- 
{SJ ' ear perturbation theory is sufficient to study the power spectrum of the primordial perturbations generated during 



> 



in 
o 



inflation, in particular to calculate the spectral index and its scale dependence. 

One way to glean more information from the CMB is to go beyond first order perturbation theory, which allows 
the study of higher order statistics, such as the bispectrum, or to calculate the amount of non-gaussianity expected 
from ones favourite early universe model. 

On super-horizon scales, i.e. scales much larger than the particle horizon, there are mainly two approaches to 
study higher order effects such as non-gaussianity: the first uses second order perturbation theory following Bardeen 
[E IS IM tm Q H, [13, EE E3, E3, EB. Ea EH H3, H3, , the second uses nonlinear theory and different 

manifestations of the separate universe approximation, either employing a gradient expansion, as originally used by 
1 Salopek and Bond |'2.1 124 . l25l l2fi| or using the AN formalism |27l l28l 123 ] 1 . In the formalism used by Salopek and 
Bond |2.'"{| the metric and the matter variables are not expanded into a power series, instead an expansion i n sp atial 
1 gradients is used. The AN formalism was introduced at linear order by Sasaki and Stewart in (3l| (see also |32ll33|) 
and relates the comoving curvature perturbation to the perturbation in the number of e-foldings and has recently 
been extended to the non- linear case |28i l2flj . In the Bardeen approach the metric and the matter fields are expanded 
in a power series in a small parameter. There has been increased interest in second order perturbation theory due to 
the prospect of new high precision data and following the papers by Acquaviva et al. @ and Maldacena Q focusing 
i-S^ ■ on the study of non-gaussianity from inflation. 

In this article we concentrate on the scalar field dynamics at second order on large scales in a universe dominated 
Q by multiple scalar fields, including metric perturbations. In linear theory the Klein-Gordon equation, including metric 
perturbations, can be written as a system of coupled evolution equations in terms of the linear field fluctuations in 
the flat gauge, the Sasaki-Mukhanov variables |2J, We show here that the Klein-Gordon equation at second 
order on large scales in the multiple field case can also be written solely in terms of the Sasaki-Mukhanov variables, 
albeit at second and first order. The Klein-Gordon equation is linear in the second order perturbations, with source 
terms quadratic in the first order perturbations (the same holds for the field equations). We also give an expression 
for the curvature perturbation on uniform density slices at second order, £2, in terms of the field fluctuations on flat 
slices at first and second order. After calculating the evolution of the field fluctuations at first and second order we 
therefore immediately get the evolution of the curvature perturbation This is similar to the AN formalism, where 
the evolution of the comoving curvature perturbation is given by the scalar field dynamics in the background, but 
which necessitates the slow roll approximation at higher order. 

We consider scalar perturbations including and up to second order in the metric and in the scalar fields using the 
Bardeen approach. We focus on scales larger than the horizon, which allows us to neglect gradient terms. We do not 
include first order vector and tensor perturbations since they do not couple to the field fluctuations on large scales; 
also scalar fields do not support vector modes and the tensor contribution from inflation is small and constant or 
decaying. The large scale focus also makes the equations more transparent and displays the relevant physics more 



1 Recently there has also been a covariant study 30]. 



succinctly, a big bonus at second order. We leave the inclusion of small scales, and vector and tensor perturbations 
to a future publication 36] . 

The paper is organised as follows: in the next section we give the field equations and the Klein-Gordon equation 
up to second order in gauge dependent form. In Section II I II after reviewing the behaviour of perturbations at first 
and second order under gauge transformations, we construct physically transparent gauge-invariant combinations. In 
Section we show how the gauge-invariant variables defined in different gauges are related to each other. We also 
discuss the construction of total variables in the presence of several fields or fluids. We review the derivation of the 
Klein-Gordon equation at first order in Section In Section lYl Al we give the field equations at second order in the 
flat gauge and in Section fVl Bl we finally derive the Klein-Gordon equation at second order. We give the second order 
Klein-Gordon equation in the single field case in Section fVIII In Section fVIIII we apply the formalism to a simple two 
field inflation model. We conclude in Section 17x1 

Throughout this paper we assume a flat Friedmann-Robertson- Walker (FRW) background spacetime and work in 
conformal time, r\. Derivatives with respect to conformal time are denoted by a dash. Greek indices, /j, v, A, run 
from 0, ... 3, while lower case Latin indices, i,j, k, run from 1, ... 3. Upper case Latin indices, /, J, K, denote different 
scalar fields. 



II. GOVERNING EQUATIONS 

The covariant Einstein equations are given by 

Gf+v = 8ttG Tpv , (2.1) 

where G^ v is the Einstein tensor, the total energy-momentum tensor, and G Newton's constant. Through the 
Bianchi identities, the field equations 1)2. l|l give the local conservation of the total energy and momentum, 

V A1 T^ = 0, (2.2) 

where V M is the covariant derivative. The energy momentum tensor for N scalar fields minimally coupled to gravity 
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<PK,n<PK,v - -j9y,v9 at 1 <PK,a<fK,0 



K=l 



g^U(<Pi, ■ ■ ■ ,¥>n) , (2.3) 



where ipx is the Kth scalar field and U the scalar field potential and (fiK.^ = §ffr- 
We split tensorial quantities into a background value and perturbations according to 

T = T + 5Ti + ^ST 2 + . . . , (2.4) 

where the background part is a time dependent quantity only To = To (77), whereas the perturbations depend on 
time and space coordinates x M = [rj,x l ], that is <5T„ = 6T n (x^). The order of the perturbation is indicated by a 
subscript, e.g. 6T± = 0(e), 

The metric tensor up to second order, including only scalar perturbations, is 

.goo = -a 2 (1 + 20! + fo) , (2.5) 

got = a 2 [Bi + ^B^j , (2.6) 

9ij = a 2 [(1 - 2V>i - ^2) % + 2E ltij + E 2 ,ij] , (2.7) 

where a = a(rf) is the scale factor, 77 conformal time, Sij is the flat background metric, <f>\ and 4>2 the lapse function, and 
ipi and ip2 the curvature perturbations at first and second order; B\ and Bi and ~E\ and E2 are scalar perturbations 
describing the shear. The contravariant form of the metric tensor is given in Appendix 151 



A. Einstein tensor 



We expand the Einstein tensor in a power series according to Eq. I|2.4[) up to second order 



= G^ v + SG^ V + 2 SG (2) V ■ 



Then the components of the Einstein tensor at zeroth order are 
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at first order 
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and at second order 
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where 0(k n ) denotes terms of order n or higher in gradients, since we are mainly interested in the large results, i.e. 
the case k — > 0. 

B. Energy- momentum tensor 

We split the scalar fields ipi into a background and perturbations up to and including second order according to 
Eq. J23}, 



The potential U = U(ipi) can be split similarly according to 



where 



U{ipi) = U + SUi + -SU 2 , 



SUi = ^ U,tp K 5wKi , 
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8tpi K 5ip 1L + ^2 U,ip K S(p2K ■ 



(2.16) 
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and we use the shorthand U <tpK = 



The energy-momentum tensor for N scalar fields with potential U(<pi) is then split into background, first, and 
second order perturbations, using Eq. 1)2. 4[l . as 



and we get for the components, from Eq. (|2.3() . at zeroth order 
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at first order 



- (e i^- 2 + tfo) , zfo - (^E ^k 2 - ^ * 

ff ( O l)0 = — \ E (^OAT^l-fr' ~ ^OA^l) - ) 

a R 

~3E(^ 1Jf -i) ' 



A" 



«?1 



E ("^oa^ia' - Voir <M ~ a 2 5U 1 



A 



5T k 



A 



ST, 



(2)j 



. A 



a 2 <W 2 



(2.20) 
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and at second order in the perturbations 

£T ( 2)o = E (^>a"W - Vtuc^^iJc' - ^oa 2 02 + Vo***? + ^ 1/f ' 2 + a 2 8U 2 ) + 0(k 2 ) , (2.25) 
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(2.26) 



6% + 0(k 2 ). (2.27) 



The — component of the energy momentum tensor gives the corresponding energy density at zeroth, first and 
second order for a universe filled by N scalar fields [jj, 



T (°o)o = -Pa > 5T n)o = ~ s Pi > ST (2)o = ~ S P2 ■ 



(2.28) 



C. Evolution of the background 

The evolution of the background is governed by the Friedmann equation given from the — component of the 
Einstein equations 



H 



8ttG / 1 
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(2.29) 



where 7i = —. The scalar fields are governed by the Klein-Gordon equation, which is from the energy conservation 
equation, at zeroth order for the ith field given by 



ip^ I + 2Hip' 0I + a 2 U^ I =0. 



(2.30) 



D. Perturbed Klein-Gordon equation 

Using the energy momentum tensor given above in Section III Bl and the energy conservation equation (|2.2|l we get 
the Klein-Gordon equation in the multiple filed case for the field I at first order 



5<pu" + 2«W + 2a 2 U, VI <h. ~ Mi^'i ~ Voi<f>'i + « 2 E u m<PK&<PiK + 0(k 2 ) = . 



(2.31) 



A 



We get the Klein-Gordon equation at second order from Eq. 12. 2|) and using Eq. i|2.31[) 

S(p 2 i" + 2H5ip 2I ' + 2a 2 U tVl (j)2 ~ 6?/>i<Vi/ - 2<j)'i5ipu + ^j^i^'i + 40ia 2 ^ U yVlipK Sip 1K 

K 

- 12^0/01^4 - 3p' 0I ip' 2 ~ ^>o/02 + a 2 X! u ,vivk 5i P2K + a 2 X] U ,<pi<pk9l S( PikS(Pil + 0(k 2 ) = . (2.32) 

if K.L 

Note, that Eqs. lEHTT]l and lf2~3"2f are in an arbitrary gauge. 

III. VARIABLES 

In this section we first calculate the changes coordinate transformations induce in the metric and matter pertur- 
bations at first and second order. We then construct physically meaningful gauge-invariant combinations from these 
variables. 

A. Gauge transformations 

We now briefly review how tensorial quantities change under coordinate transformations 0, . We consider two 
coordinate systems, and a; M , which are related by the coordinate transformation 

a ^ = e? A afx^ j (3.1) 

where £ A is the vector field generating the transformation and £ M = £j + |£ 2 + 0(e 3 ). Equation (|3.1(l can then be 
expanded up to second-order as 

^ = ^ + tf + \{tf,Ai+$) ■ (3-2) 
A tensor T transforms under the change of coordinate system defined above as 

T = e £ « A T, (3.3) 



where denotes the Lie derivative with respect to £ A . Splitting the tensorT according to Eq. (|2.4() we get 

5Ti - 5Ti+L €l T 

5T 2 = 5T 2 + L e2 T +L| i T +2L ei OT 1 . (3.4) 

Under a first-order transformation ^" = (ai,[3 1 l ), a four scalar such as the energy density, p, or the scalar field tp, 
transforms therefore as 

dpi = Spi + p' ai , (3.5) 
and or the first order metric perturbations we have 

■01 = ipi - Hai , (3.6) 
0i = 0i + Hai + a[ . (3.7) 
At second order, writing — (a 2 ,/3 2 l ), we find from Eq. I|3.4(l that a four scalar transforms as 

Sp 2 = 5p 2 + p' a 2 + ai [p'o a i + Po a i + 28 Pi] 

+ (2S P i+p' a 1 )^(3 i ;, (3.8) 

and the second order metric perturbations transform as 

4>2 = V^2 - Ha 2 - on [Hon' + {W + 2H 2 ) a t - 2ip[ - AU^i] 

-(Wai-20i) fe /3 1: fc , (3.9) 
4> 2 = 02 + Uoi 2 + a 2 ' + ai [ai" + 5Hai + (W + 2H 2 ) on + AHfa + 2(p[] 

+ai' (2a/ + 40i) + (3i, k {on' + Hai + 2fa) k + f3[ k («i - 2B X - 0[) k . (3.10) 



From Eqs. (|3.8[1 . I|3.9|l . and l|3.10|l we see that on super horizon scales, where gradient terms can be neglected, the 
definition of the second order perturbations in the "new" coordinate is independent of the spatial coordinate choice 
(the "threading") at second order in the gradients. It is therefore sufficient on lar ge s cales (at 0{k 2 )) to specify the 
time slicing by prescribing a\ and 0:2, in order to define gauge- invariant variables [lj,|2g. The procedure to neglect 
the gradient terms, is explained in detail in Ref. |28[ . For the approximation to hold one assumes that each quantity 
can be treated as smooth on some sufficiently large scale. Formally one multiplies each spatial gradient di by a 
fictitious parameter fc, and expands the exact equations as a power series in k, keeping only the zero- and first-order 
terms, finally setting k = 1. 

Note, in the following sections we shall omit the symbol 0(k n ) denoting the order of the gradient terms neglected, 
assuming that if not stated otherwise the equations are valid to 0(k 2 ). 



B. Gauge-invariant combinations 

Using the transformation behaviour of the perturbations derived in the last section we can now construct com- 
binations of these variables that do not change under gauge transformations, i.e. gauge-invariant variables. The 
combinations given below have definite physical meaning, such as the curvature perturbation on a hypersurface of 
uniform density. 

We define the various time slicings and then substitute the results into the transformation equations to arrive at 
the gauge-invariant variables. Whereas the first order results are valid on all scales, we only consider the large scale 
case at second order. 



1. First order 

Hypersurfaces of uniform /-field, i.e. Sifiu = 0, are given from Eq. (|3.5|) by the temporal slicing 

a\ = — . (3.11) 

<Poi 

The curvature perturbation on uniform field hypersurfaces |37t l3^ | , or comoving curvature perturbation, can then be 
defined for each field using Eq. (|3.6|) as 

n 1I = ip 1 + n^r L . (3.12) 

</>0J 

Flat slices are defined as ipi = and we therefore get from Eq. H3.6JI the first order time shift 



t 

The Sasaki-Mukhanov variable [34l I35I ]. or the field fluctuation on uniform curvature hypersurfaces, is then given by 



«i = fi- (3-13) 



Qu = S VlI = Sifiu + ■ (3.14) 

rt 

The density perturbation on uniform curvature hypersurfaces is defined as 

<W = S Pla + Vi . (3-15) 

where po a and 8p\ a are the energy density of the a-fluid in the background and at first order 2 . The curvature 
perturbation on uniform a-fluid density hypersurfaces is at first order defined as 

Ci a = -A-n 6 -^. (3.16) 



The lapse function on flat slices is from the definition l|3.7|l and Eq. I|3.13|l given as 



^i+fl-SjUi. (3-17) 



2 Greek indices from the beginning of the alphabet, a, /3, 7 will be used to denote different fluids. 



2. Second order 



At second order uniform field slices are, from Eq. (|3.8|) . defined by the temporal shift 

Ot-2 = 7- ( 8<f2i —Sifn'Sipu) , 

<Poi V Vol ) 



(3.18) 



on large scales. The curvature perturbation on uniform /-field slices is therefore on large scales at second order, from 
Eq. ^ and using (|3~TT|) and fHM 14] 



I 2 

<Foj 



Voi 



(3.19) 



The time shift that defines the fiat slicing is at second order, from Eq. (|3.9(l . on large scales given by 



n 



n 



(3.20) 



Then the field fluctuation on uniform curvature hypersurfaces, or Sasaki-Mukhanov variable, at second order for the 
Jth field, from Eq. and definitions of the time shifts lf3~T3|) and IpHO")) is 



<Poi .,. , ( w 



0,21 = Slfi2I = Slfi2I + -j^1p2 + ( 77 



2H<p' 0I + tp' \ - j£<Poi 



(3.21) 



The second order lapse function in the flat gauge is given from Eqs. H3.10(l . Ij3.13|l . and l|3.2()(l as 

nil \ i / T-/' \ / T-l 1 T-/" T-t 1 ^ N 



(3.22) 



The curvature perturbation on uniform a-density hypersurfaces is at second order on large scales from Eq. i|3.9|) 
and using l|3.11|l and (|3.18|) given by [13] 



(2a = ^ + ^-8p 2a - 2^L.5p[Jp la ~ 2^ M + 2H^) + 6 -^4 (n^ -H'- 2H 2 

POa Poa POa Po a \ Poa 



(3.23) 



and from Eq. (|3.8|) and the definitions of the time shifts (|3.13|) and 13.20fl we get the second order density perturbation 
on uniform curvature hypersurfaces [l4| 



(3.24) 



IV. RELATING VARIABLES 

In this section we relate the gauge-invariant variables defined above on different hypersurfaces to each other. Using 
the field equations we then define total perturbations from the individual field and fluid ones, where possible. 

A. First order 

From Eqs. (|4.5|l and (|3.14(l the Sasaki-Mukhanov variable of the field / is related to the curvature perturbation on 
uniform field hypersurfaces at linear order simply by 



(4.1) 



From Eq. I|3.15[l and i|3.16[l we find that the curvature perturbation on uniform density hypersurfaces is related to 
the density perturbation on uniform curvature hypersurfaces simply as 

ClQ = -H%. (4.2) 

The total density perturbation at first order in the flat gauge is given in terms of the density perturbations of 
individual fluids simply by 

Spi = Spi a , (4.3) 

a 

which allows us to relate the total curvature perturbation to the individual fluid curvature perturbations by 

Ci = £ ^Ci« • (4.4) 

We can define the total comoving curvature perturbation (i.e. relative to the average fluid velocity or an "average 
field" ) using the — i component of the field equations as 

l^K VoK I 

The total curvature perturbation on uniform density hypersurfaces in terms of the field fluctuations on flat slices is 
given by 

Ci = ~ H , a X>o*r&Ti- (4-6) 
L-ih Pol k 

B. Second order 

The Sasaki-Mukhanov variable at second order defined above in Eq. (|3.21|l can be expressed in terms of the curvature 
perturbations on uniform field hypersurfaces, IZu and 7£.2/, 



2 



Z^jKiKu ■ (4.7) 



Similarly we can express the curvature perturbations on uniform field hypersurfaces in terms of the Sasaki-Mukhanov 
variables at first and second order and get, 

U 2I = ^-Q„- + % (n^ -H>- 2H 2 ) . (4.8) 

To define the total comoving curvature perturbation at second order 72-2^ in terms of the 'R.ii we would need the 
— i Einstein equation at this order in an appropriate form, that is without gradients. Since it is not clear how to 
arrive at this form of the — i Einstein equation (without imposing slow roll, see Section IVl Al below ) we shall leave 
the definition of 7Z21 open for the moment and shall return to this issue in a future publication [^a]. We note however, 
that the definition of the total comoving curvature perturbation at second order 1Z2 V is not a problem in itself, since 
it was already shown in Ref. that on large scales TZ2 V and (2 coincide. The definition of the total curvature 
perturbation (2 in terms of the field fluctuations is given below in Eq. (|4.16(l . and we can therefore get 72.2^ from £2 
if we should need it. 

The definition of "total" quantities from quantities defined for a specific field or fluid is not more problematic 
at second order than at first order, it is a mere question of having the "right" equations. As an example we shall 
now digress slightly from the main theme of this paper, multiple scalar fields, and analyse the definition of the 
total curvature perturbation at second order in the case of multiple fluids, which is unlike the definition of the total 
comoving curvature perturbation, straight forward 3 . 



3 The rela tion between scalar fields and fluids at first order, in particular how fields can be treated as fluids, has been studied in detail in 
Ref. [33. 



We can relate the curvature perturbation on uniform a-fluid hypersurfaces, Eq. (|3.23() . to the density perturbation 
on flat slices, Eq. IpHI)) . 



U ~ W ~ / ~ <5pi Q 

-C2q = -^—0p 2a - Z—^Opiu dp la H ^— 2" 
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and similarly the density perturbation of the a-fluid at second order on flat slices in terms of Cia and (2c 
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The total density perturbation at second order on flat slices is given in terms of the density perturbations of 
individual fluids on large scales simply by 



5p2 = ^2 3P 2 °> ' 



(4.11) 



which allows us to write the total curvature perturbation at second order to the individual fluid curvature perturba- 
tions, using Eq. I|4.1U[) . as 



c 2 = E% 2Q -(i+3c^)c 1 2 +E 
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where we used the background energy conservation equation for the a- fluid, Eq. QA7J) . given in Appendix c 2 a is the 
adiabatic sound speed of the a-fluid, and the total adiabatic sound speed is denoted c£. 4 
The relative entropy (or isocurvature) perturbation at first order is defined as |4ll |42( 

Si a p = 3 (Ci - Ci/s) , (4-13) 
which can be used to rewrite Eq. Ij4.12|l in terms of adiabatic and non-adiabatic quadratic first order terms as 
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(4.14) 



The evolution equation for the relative entropy perturbation at first order, <Si Q7 , was given in |39l |. We see that there is 
a purely non-adiabatic contribution to £2, quadratic in <Si Q7 . Note that for purely adiabatic first order perturbations 
Eq. Ij4.14jl simplifies to 



f _ \ " Poa /■ 
„ ^0 



(4.15) 



Turning back to scalar fields, the curvature perturbation (2 in terms of field fluctuations on flat slices at first and 
second order is, using Eq. I|3.23[) evaluated for a single fluid, the energy densities in terms of scalar fields Eq. 
and the definitions for the Sasaki-Mukhanov variables Eqs. I|3.14fl and l|3.21(l . given by 
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(4.16) 



4 The rela tion between the total £2 and the curvature perturbations of the individual fluids, C,2a, has been given in the two field case in 
Ref. E3. 



where 8U2 is given from Eq. Ij2.19|) above evaluated on fiat slices, 



Qil+J2 U ^kQ2K, (4.17) 

K.L K 

and the adiabatic sound speed for a universe filled by N scalar fields is given by 

and po is the background energy density defined in Eq. I|2.28[) above. 

Equation l|4.16[) can be readily evaluated either analytically or numerically given the Sasaki-Mukhanov variables at 
first and second order, Qu and Qii- However, in Section IVf f I below we use — i Einstein equation at second order 
in the single field case and and in Section IVf f f I we use the slow roll approximation in a particular model to simplify 
Eq. i|4.16[l further, without the time derivatives of the Sasaki-Mukhanov variables. 

V. GOVERNING EQUATIONS IN THE UNIFORM CURVATURE GAUGE AT FIRST ORDER 

We begin this section by giving the field equations at first order in the uniform curvature gauge and then review 
the derivation of the Klein-Gordon equation for multiple scalar fields at first order on large scales. 

A. Field equations 

The components of the Einstein tensor and the energy-momentum tensor are given in Sections [O and ITlBl 
Substitution into Eq. I|2.1|l then gives for the — component of the Einstein equations gives at first order 

2a 2 U ^i+Y, i PaKQiK' + a 2 6lh =0, (5.1) 

K 

and from the i — j component of the Einstein equation we get 



a ~i 
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i -87rG^2ip' 0K Q 1K ' = 0. (5.2) 



K 



At first order the — i components of the Einstein tensor, Eq. l|2.11|l . and energy momentum tensor, Eq. l|2.23|l . 
are linear in the spatial gradients (by definition the background quantities are only time dependent), which allows us 
to write the — i Einstein equation without gradients as 



-<£i-47rG^y 0K g lif = 0. (5.3) 



K 



Combining Eqs. JO} and l|5.3|l we get 

E ^kQik = -a 2 E (u VK + ^EWojt) Qik , (5.4) 
relating the time derivative of the field fluctuations to the field fluctuations themselves. 

B. Klein Gordon equation 

At first order the Klein-Gordon equation for the field /, on flat slices, is from Eq. I|2.31|l given by 

Qu" + VHQ11' + 2a 2 U Vl fa - p' 0I fa + a 2 E U mvx QlK = , (5.5) 

K 



which can be simplified using Eqs. (|5.1() and (|5.2() . to give 



Qu" - WQu' • • ( ^ + f'oi^- ) E ^'okQik + a 2 E ( U mVK - -f^f, ,(;,,; j =0. (5.0) 
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Equation 1)5.6)1 can be further rewritten using Eq. (|5.4|) above, to give 0, Q 01 
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Cu" + 2HQ U ' + E 



A' 



87rG / 2 , / ¥>0A'\ V 



QlA = . 



(5.7) 



which now displays the "canonical" time derivatives (in conformal time) d 2 /drj 2 + 2Hd/dri, but has a rather involved 
mass term. 

VI. GOVERNING EQUATIONS IN THE UNIFORM CURVATURE GAUGE AT SECOND ORDER 

In this section we first give the field equation in the uniform curvature gauge at second order, highlighting problems 
arising from the — i equation at second order and possible ways to solve them. We then give the Klein-Gordon 
equation in the multiple field case on large scales in terms of the Sasaki-Mukhanov variables. 



A. Field equations 

In this section we present the field equations in the uniform curvature gauge at second order on large scales. The 
components of the Einstein tensor and the energy-momentum tensor are given in Sections III Al and III Bl Substitution 
into Eq. 1)2.1)) then gives for the — field equation in the fiat gauge at second order is 



2a 2 C/ o (40i 2 -; 

and, using Eq. 1)5.3)1 and Eq. (|6.1|) above, the i — j Einstein equation can be rewritten as 



E (<PokQ2K - ^OkMk + Q'lK 
K 



a 2 5U 2 



(6.1) 
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= 8ttG 
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The — i Einstein equation at second order is given by 

H - 40i?i, 4 ) - 4ttGE (v'ok QiK,i + ZQik'Qika) + 0(k 3 ) = , 



(6.2) 



(6.3) 



A" 



where we used Eq. 1)5.3(1 . The left hand side of Eq. 1)6. 3JI can be rewritten as TL (j)2 — 20i ) , with the spatial 

derivative outside the brackets. However the right hand side of Eq. 1)6.3)1 can't be written, in the general case, as an 
overall gradient, as can be seen above: whereas the Qik,i term is multiplied by a background quantity, similar to the 
first order case, the QiA'QiA.i term can't be written as an overall gradient and the — i Einstein equation at second 
order cannot be brought into scalar form immediately. 

In order to recast Eq. 1)6.3)1 in a more useful form, without the gradients, we have several possibilities: 

• The easiest solution is to require Q\k = const or Q\k = 0. Unfortunately this case is not particularly 
interesting, since one of the most interesting applications of second order theory is the study of non-gaussianity, 
which necessitates the inclusion and evolution of the first order and in particular the first order squared terms. 
We shall therefore not pursue this option further. 



• The next case is to assume we only have one field to consider and we shall study this case in detail in Section 
IVHI Actually, we need only one field at first order, but to assume a single field at first order and multiple field 
perturbations at second order seems to be rather contrived. 



• Another possibility is to use the slow roll approximation such that Eq. (|5.7(l can be rewritten as Q' u oc f(ipoj)Qu, 
where f(<fioj) is a function of the background fields, which allows as in the single field case to replace Q' u Qn 
by 6? T in Eq. (|6.3|1 (for a recent detailed exposition of the slow roll approximation in the multi-field case see 
e.g. 22]). We shall illustrate this particular case, i.e. using the slow roll approximation, by studying a simple 
two-field inflation model in Section lYlIII 

• Finally, we can take the divergence of Eq. Q6.3|l . which we shall do next. 
Following Ref. @ we take the divergence of Eq. (|6.3() and get 

nv 2 J 2 - mUiV^ + ^l^ 



- 4ttG^ [p' k^ 2 Q2K + 2Q' 1K V 2 Q 1K + 2Q' 1Kj Q 1K j ] + 0(k 3 ) = 

K 

Using now the inverse Laplacian, defined as V~ 2 V 2 X = X, Eq. 16.411 can be rewritten as 

H (h - Hh) - ^GJ^HkQsk + 2Q' 1K Q 1K ] - R(Ti, Si/) + 0{k 3 ) = , 

K 



where we define R{(f>i, Qu) as 

R(Tu Qu) = 4W 



v- 2 4>iV 2 4>i 



n,k<Pi, 



8nGY, [v- 2 Q' 1k V 2 Q 1 k + V- 2 (o! 1K<j Q lK ; 



(6.4) 



(6.5) 



(6.6) 



But since neither the effect of R((f>i, Qu) on the evolution of the field fluctuations nor its large scale limit is clear, we 
shall not use this form of the — i Einstein equation below. 

Not being able to make use of the — i Einstein equation is inconvenient, but since it is only a constraint equation 
and therefore redundant, we are able to get a closed system of equations without using it, as shown in Section IVI Bl 
below, which is sufficient to get the evolution of £2 from Eq. (|4.16|l . However, the — i Einstein equation does allow 
us to rewrite the Klein Gordon equation in a more compact form with canonical time derivatives (see above the first 
order case, Section and below the single field second order case, Section IVlTll . 



B. Klein Gordon equation 

In this section we derive the Klein-Gordon equation in the flat gauge for N scalar fields on large scales. 

At second order the Klein Gordon equations on flat slices is for the field /, from Eq. (|2.32|l , on large scales given by 

Q21" + 2HQ2I 1 + 2a 2 U^J 2 ~2j l Q 1I ' + A v '^ x Ji +Aj ia 2 Y,U^ lVK QiK 

K 

- <Pai4>2 + a 2 ^ U,ipi VK Q2K + fl2 ^ U^jpkvl QikQil = ■ 



(6.7) 



A' 



K,L 



We can rewrite Eq. H6.7(l above, using the field equations at first and second order given above in Sections IV Al and 
IVI Al to substitute for the lapse functions at first and second order, cf>i and cf>2, and get 



Q21" + 2HQ2i' + a 2 J2 




U :Vl U, i( p Kl p L + 



Q2K 
16?rG 



H 



f + ^0i)E(^kQ2k'+Q 1 k' 2 ) 

^2 {^wpk - ^ U ,<pi U ,<pk^ ^ol| QilQik 



<Pok) QikQu = 0. 



(6.8) 



This is the gauge-invariant Klein-Gordon equation in the uniform curvature gauge for iV minimally coupled scalar 
fields on large scales in terms of the field fluctuations in the flat gauge at first and second order, Qu and Q,2i- It is 
linear in the second order variables, but quadratic in the first order ones. 



VII. KLEIN GORDON EQUATION AND CURVATURE PERTURBATION: SINGLE FIELD CASE 



Having dealt with the multi-field case in the previous section, we now turn to the single field case. The general 
form of the Klein Gordon equation at second order, Eq. (|6.8|) . reduces in this case to 
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V'o) Qi'Qk = 0. 



(7.1) 



In the single field case we can use Eq. H5.4(l to rewrite the — i Einstein equation at second order, Eq. (|6.3|) . as 

'2^ 



H (fa - 2 fa ) 



4ttG 



( (p H 



+ O{k 6 ) = 0, 



(7.2) 



which allows us to immediately get rid of the spatial gradient in a similar fashion as in the first order case. From 
Eqs. (f f> . ip and Q7.2J1 we then get the useful relation, similar to Eq. (|5.4(l at first order, 
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(7.3) 



Substituting Eq. (|7.3(l into Eq. (|T. 1|> we arrive at the single field Klein Gordon equation at second order, 



Q 2 " + 2HQ 2 ' + a 2 
+ a 2 {U. 



U., 



16ttG 
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(7.4) 



now with the time derivatives in canonical form, which is nice. 

In the single field case the curvature perturbation on uniform density hypersurfaces at first order in terms of the 
field fluctuation on the flat slice is given, from Eq. Q4.6[l . as 



Ci 



4s, 

Pa 



and at second order we get from Eq. (|4.16() , using Eq. 17. 3|) 
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(7.5) 



(7.6) 



Using Eqs. (|7.4|) and 117. 611 s hould be sufficient to show the conservation of C2 on large scales in the single field case, 
which was shown in Refs. |l8l |23| | , and we will return to this issue in a future work |36| . 



VIII. A SIMPLE APPLICATION USING SLOW ROLL 

We can now apply the formalism derived in the previous sections to the simple two field model of Ref. which 
has subsequently been studied in Refs. [2fl , I29I ] and most recently in Ref. (f° r earlier work on two-field inflation at 
linear order see e.g. ^3)- Using the slow-roll approximation we calculate the curvature perturbation at second order, 
£2, in terms of the field fluctuations. 



The potential is given as 



U = U + -mlvf + -m%<p\ , 



(8.1) 



with the first term, Uo — const, dominating and where mj denotes the mass of the 7th field. Furthermore we assume 
slow roll and 



<A)2 = , ^02 = • 

The background Friedmann constraint, Eq. I|2.29[) . then simplifies to 
and the background Klein-Gordon equation, Eq. Il2.30|) . gives 

2 2 

Vol + ~^f^ 01 = °- 



•2) 



(8.3) 



(8.4) 



From Eq. (|4.6|l we immediately get the curvature perturbation at first order in terms of the field fluctuations on 
flat slices in this model, 



Ci = 8ttG- 



Uq 



-Qi 



•5) 



The calculation of the second order curvature perturbation in terms of the field fluctuations is slightly more involved. 
As pointed out in Section llV Bl in order to get ( 2 solely in terms of the field fluctuations without their time derivatives 
we have to use the — i Einstein equation at second order, and, as discussed in Section IVf Al to get the latter in a 
useful form we need the first order Klein-Gordon equations. 

The perturbed Klein-Gordon equation at first order on flat slices, Eq. I|5.5|l . reduces in the slow roll limit for the 
model specified by Eqs. I|8.f Jl and (|8.2|l to 
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1 Sn = o, 
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where Qn and Q12 are the field fluctuations in the flat gauge, the Sasaki-Mukhanov variables, at first order for the 
two fields. 

The — i equation at second order on flat slices, Eq. I|6.3J) . simplifies then using Eqs. Ij8.6|l to 



U (02 - 2$) = 4ttG 



Vox 0,2 
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Rewriting Eq. (|6.1|l . the — Einstein equation at second order, in a similar fashion we arrive at the useful relation 
J2 (wkQ 2 k + Q'ik) + a 2 SU 2 = a 2 {m\ Vol Q 21 + mjQ 2 , + m 2 Q 2 2 ) , (8.8; 



K 



which upon substitution in Eq. I|4.16|l yields 
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(8.9) 



the curvature perturbation on uniform density hypersurfaces at second order in terms of the field fluctuations on flat 
slices for the model specified above. 

This expression for ( 2 , Eq. I|8.9|l above, agrees with the one found by Lyth and Rodriguez in Ref. [23] using the 
AiV formalism, if we take into account that £2 = C2LR + 2£i 2 [2^]. In particular we find that (2 does not contain any 
non-local terms. 

However, the curvature perturbation £2 found here disagrees with the expression for the second order curvature 
perturbation found by Enqvist and Vaihkonen in Ref. [jj]], although the order of magnitude estimate published 
subsequently by Vaihkonen |46| seems to agree with the result found here. 



IX. DISCUSSION AND CONCLUSIONS 



In this paper we have presented the Klein-Gordon equation for multiple minimally coupled scalar fields at second 
order in the perturbations in a perturbed FRW background on large scales. We have shown that using suitable gauge- 
invariant variables, namely the field fluctuations in the flat gauge or Sasaki-Mukhanov variables, the Klein-Gordon 
equation at second order, Eq. I|6.8[l . can be written solely in terms of these variables, as in the first order case, and is 
linear in the second order variables, but has source terms quadratic in the first order field fluctuations. 

We have also given the relation between gauge-invariant quantities in different gauges, and hence on different 
hypersurfaces, which at second order is non-trivial. In particular we give the curvature perturbation on uniform 
density hypersurfaces, £2, in terms of the Sasaki-Mukhanov variables of the individual fields, Qu and Q21 in Eq. i|4.f 6fl . 
We calculated (2 for a particular two-field model durin g slo w-roll inflation, Eq. 1)8. 9)) . and found excellent agreement 
with the expression derived using the AN formalism |29j |. In particular we find using second order cosmological 
perturbation theory that there are no non-local terms in the expression for (2 in this model. 

Having an expression for £2 solely in terms of the field fluctuations allows us to get the evolution of £2 directly 
from solving the Klein-Gordon equations at zeroth, first, and second order, Eqs. 1)2.30)1 . 1)5. 7|l . and (|f> .8|> . respectively, 
together with the Friedmann constraint, Eq. 1)2.29)1 . Alternatively we could have calculated the non-adiabatic pressure 
due to the multiple scalar fields first and then solve an evolution equation for £2 [iJ, l39i |41| , but using Eq. (|4.f 6|) 
eliminates the integration of the £2 evolution equation and is therefore simpler, since the Klein-Gordon equations have 
to be solved in both cases. 

Given suitable initial conditions for the Klein-Gordon equations governing the dynamics of the fields at zeroth, 
first, and second order the equations can be integrated numerically. The equations can be solved order by order, the 
background and first order fields acting as source terms for the second order equations, since we do not consider back- 
reaction of the higher order perturbations on the lower order ones. Solving the second order equations numerically 
will therefore be very similar to solving the first order system. 

Alternatively, employing a slow roll approximation at all orders, the equations can be approximated analytically. It 
will be particularly interesting to compare the evolution of £2 from a slow roll version of the Klein-Gordon equation, 
Eq. 1)6.8)1 , to the AiV-formalism [2£j , which also uses slow roll and also gives the evolution of the curvature perturbation 
directly from the solution of the Klein-Gordon equations, without having to integrate an evolution equation for £2- 

The next step will be the extension of the formalism presented here to small scales, i.e. deriving a formalism valid 
on all scales. This will also be an advantage compared to other approaches based on the separate universe paradigm, 
such as the AiV-formalism, which can not be extended to small scales since it is by assumption only valid on scales 
of order or larger than the horizon. First steps towards extending the formalism are under way at present, based on 
the work presented here |3^ . 
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APPENDIX A: BACKGROUND FIELD EQUATIONS 



The Friedmann equation is given from the — component of the Einstein equations 



8ttG 



(Al) 
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a po , 



where H=^. 



From the i — j component we find 




(A2) 



The two previous equations can be rewritten as 



—a 2 (po + 3P ) , 
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or alternatively as 
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Also useful are 

2^V _^ =47 rGa 2 (P + Po) , 
\ a / a 

where the right hand side in the case of scalar fields simplifies to a 2 (Pq + po) = J2k Voif 2 j an d 

W + 2H 2 = AnGa 2 (p -Fit), 
where the right hand side in the case of scalar fields simplifies to (po — -Po) = 2£/o- 
The background energy conservation for the a-fluid is given by 0, EE1 

P[)a = ^ 37Y {poa + Poa) + Qa , 
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where po ai Po a , and Q a are the energy density, the pressure and the energy transfer to the a- fluid. Note that the 
energy transfer defined here is related to the one define in 39] , Q a , by aQ a — Q a . The adiabatic sound speed of the 
a-fluid is 
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related to the total adiabatic sound speed by 
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For more details on the multi-fluid formalism see 



APPENDIX B: THE METRIC TENSOR 



The metric tensor up to second order, including only scalar perturbations, is 

.goo = -a 2 (1 + 20i + cj> 2 ) , 
f 



and its contravariant form is 
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